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ABSTRACT 

Microlensing is a powerful tool for studying stellar atmospheres because 
as the source crosses regions of formally infinite magnification (caustics) 
the surface of the star is resolved, thereby allowing one to measure the 
radial intensity profile, both photometrically and spectroscopically. However, 
caustic crossing events are relatively rare, and monitoring them requires 
intensive application of telescope resources. It is therefore essential that the 
observational parameters needed to accurately measure the intensity profile are 
quantified. We calculate the expected errors in the recovered radial intensity 
profile as a function of the unlensed flux, source radius, spatial resolution 
the recovered intensity profile, and caustic crossing time for the two principle 
types of caustics: point-mass and binary lenses. We demonstrate that for both 
cases there exist simple scaling relations between these parameters and the 
resultant errors. We find that the error as a function of the spatial resolution 
of the recovered profile, parameterized by the number of radial bins, increases 

3/2 1/2 

as Nft , considerably faster than the naive expectation. Finally, we 

discuss the relative advantages of binary caustic-crossing events and point-lens 
events. Binary events are more common, easier to plan for, and provide more 
homogeneous information about the stellar atmosphere. However, a sub-class 
of point-mass events with low impact parameters can provide dramatically 
more information provided that they can be recognized in time to initiate 
observations. 
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1. Introduction 

Originally proposed by Paczyriski (1986) as a method to detect the presence of 
massive compact objects in the halo of our Galaxy, microlensing has increasingly been 
recognized as a tool for studying a broad range of astrophysical phenomena. Various 
applications include detection and characterization of binary and planetary systems (Mao 
k. Paczyhski 1991; Gould Sz Loeb 1992), reconstruction of the stellar mass function down 
to masses below the hydrogen burning limit (Paczyhski 1991; Griest et al. 1991; Gould 
1996), measurement of the rotation speed of giants (Gould 1997), measuring the transverse 
velocity of galaxies (Gould 1995a), and probing the central engines of quasars (Gould & 
Gaudi 1997). Recently, Valls-Gabaud (1994, 1997) and Sasselov (1996) have proposed 
using microlensing to study stellar atmospheres. Here we analyze this application in detail. 

Currently, three collaborations (MACHO, Alcock et al. 1997a; OGLE, Udalski et al. 
1997; EROS, Ansari et al. 1996) have ongoing projects that survey the Galactic bulge with 
roughly nightly sampling in order to detect microlensing events. Over 60 events per year 
are being detected. These data are being analyzed real time which has allowed MACHO 
to issue "alerts", notifications of ongoing events detected before the peak. Two follow-up 
collaborations (PLANET, Albrow et al. 1997, 1998a; GMAN, Alcock et al. 1997b) have 
formed in order to monitor these alerts around the clock with high precision and high 
temporal sampling with the aim of detecting (and further alerting on) deviations from 
the standard microlensing light curve, such as would be expected from binary lenses, 
binary sources, finite sources, and parallax. These deviations are useful in that they can 
provide additional information about the lens and/or source. In addition, a collaboration 
has formed to conduct spectroscopic monitoring of alerted events in order to study the 
source stars in detail and has observed several events to date (Lennon et al. 1997). Thus, 
the prospect for the real-time detection and monitoring of light curve anomalies both 
photometrically and spectroscopically is promising. 

All gravitational microlenses have caustics, defined as the set of points in the source 
plane where the magnification of a point source is formally infinite. When a finite source 
crosses a caustic, the gradient of the magnification over the source is large, and therefore 
different parts of the source are magnified by different amounts. Hence, the source is 
partially resolved. Different parts of the source are resolved at different times during 
the caustic crossing, and thus, by taking a series of measurements during the course of 
the crossing, one can recover the intensity profile of the source. Several workers realized 
that finite source effects could be useful for breaking or partially breaking the degeneracy 
among microlensing event parameters (Gould 1994; Nemiroff & Wickramasinghe 1994; 
Witt & Mao 1994; Maoz & Gould 1994) and that variations in the surface profile could 
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be exploited to this end (Witt 1995; Loeb & Sasselov 1995; Gould & Welch 1996). 
However, Valls-Gabaud (1994, 1997), Sasselov (1996), and Heyrovsky, Sasselov, <fe Loeb 
(1998) proposed to exploit the same effects for radically different aims. Instead of using 
finite-source effects to learn more about the lens, they sought to learn more about the 
source. The basic idea is as follows: Imagine that one could image separately different 
annular rings on the surface of a star. In effect one would be sampling different depths of 
the photosphere. Since the temperature varies as a function of depth, the broad spectral 
energy distribution would change with annulus, with more blue light near the center 
(greater depths) and more red light near the outer limb (lesser depths). That is, the star 
would be limb-darkened, and more in the blue than the red. Since different spectral lines 
form at different depths, one would expect that the detailed spectral profile would vary as 
a function of annular radius. Hence, the entire atmosphere could be studied as a function 
of depth by resolving the two-dimensional (radius and wavelength) spectral profile of the 
star. Currently, it is only possible to study stellar atmospheres in this way for the Sun and 
eclipsing binaries. However, since the surface of the source star is partially resolved during 
a microlensing caustic crossing, one can also probe the atmospheres of the source stars for 
these types of events. 

Both Valls-Gabaud (1994, 1997) and Sasselov (1996) used specific stellar atmosphere 
models to construct broad-band and spectral line brightness profiles, and then used these 
profiles to predict in detail the variations in the broad-band color or equivalent width 
of specific lines that one would expect during the course of a point-mass microlensing 
event. Sasselov (1996), and Heyrovsky et al. (1997) also consider the effects of star spots 
on the microlensing light curve. In addition, they compared their predictions to MACHO 
Alert event 95-30, a point-mass caustic crossing event for which spectra were taken 
during the course of the crossing, and for which variations in the optical TiO bands were 
detected (Alcock et al. 1997b). All of these authors predict that the color and spectral-line 
variations during the caustic crossing should be significant and note that this provides an 
entirely new method of studying stellar atmospheres. 

Although caustic crossing events are in principle useful for studying the atmospheres 
of stars, these events are rare, and they typically last for only about 7 hours (for a giant 
source). For this method to be successful, it is essential that observers have a clear sense 
of what can be accomplished with these events, since substantial telescope resources are 
likely to be expended. To this end, we approach this topic from another perspective. 
We quantify the intrinsic ability of both point-mass and binary lens microlensing caustic 
crossings to recover the radial variation of the intensity of the star for any arbitrary 
wavelength, and hence, for any spectral line. Specifically, we calculate the fractional error 
in the recovered intensity profile as a function of the unlensed flux of the source, the 
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duration of the measurements, the size of the telescope, the magnification of the event, 
and the spatial resolution of the recovered intensity profile. This information will be useful 
to observers in making rational decisions about which events to follow, and what resources 
are required to address specific questions. 



2. Formalism 



2.1. Lens Geometry, Magnification, and Caustics 



Microlenses come mostly in two flavors: point lenses and binary lenses. The 
magnification structure of these two types of lenses is essentially different, so we consider 
them separately. 

The flux of a point source being microlensed by a point lens, F, can be expressed in 
terms of five parameters, Fq, the unlensed flux of the source, B, the flux of any unresolved 
background light, to, the time of maximum magnification, (5, the impact parameter of the 
event, and t e , the time scale of the event. These are related by F = FoA(t) + B, where 
A(t) is the magnification, which is itself a function of the parameters to, t e , and (3, and is 
given by, 

1/2 

(2.1) 



A[x(t)]= ( X 2 ^, 1/2 , x(t) = 
x(x z + 4) i / / 



p 2 + (t~to 



,2 



Here x and (3 are in units of the Einstein ring radius, R e , given by, 

R e = ^GMD ol (1-Dol/Dos)] 1/2 ^ , 22) 
e c 2 

where M is the mass of the lens, and Dol and Dos are the distances to the lens and source 
respectively. The time scale is related to the Einstein ring radius by t e = R e /v, where v is 
the transverse speed of the lens relative to the observer-source line of sight. Note that for 
x< 1, A(x) ~ x^ 1 , and that the magnification therefore diverges for a point source. Thus 
the point x = corresponds to the caustic for the point-lens case. Consider a source of 
radius R*. The radius of this star, projected onto the lens plane, in units of R e is given by, 

For M = 0.3M Q , D l/D 0S = 0.75, D OS = 8kpc, and R* = 10R®, p ~ 0.02. For a source 
of this size, equation (2.1) remains valid for the majority of the event. However, when the 
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lens comes within ~ 2p of the source, i.e., for events with f3~ 2p, the magnification deviates 
from equation (2.1), and the finite size of the source must be considered. The time scale for 
this deviation is roughly the crossing time for the source, t C ross = pt e — 7(R*/WR & ) hours 
for typical bulge events. For these events, the source parameter p enters into the calculation 
of the magnification. The probability of such an event is ~ (p), where (p) is the average 
radius of the source stars being monitored. If 100 giant-source events were detected during 
a year, then the lens would transit the source for about 5 of them and these would exhibit 
finite-source deviations (Gould 1995b). 

The magnification of a point source being lensed by a binary depends on the same 
parameters as the point lens, to, t e , and (3, along with the additional parameters b, the 
separation of the lenses in units of R e , q, the mass ratio of the lenses, and 9, the angle 
between the axis of the binary and the trajectory. Unfortunately, the dependence of the 
magnification on these parameters in the general case has no analytic form. Since the 
method of determining the magnification in the general case has been described elsewhere 
(see, e.g. Witt & Mao 1995), and since we are only interested in caustic crossings, we will 
employ a simplified formalism for these crossings. For nearly equal mass binary lenses, the 
caustics are composed of curved lines (called folds) whose extent is of order the Einstein 
ring radius. Thus the probability of encountering a fold caustic during a binary-lens 
event is almost unity. For a point source, the excess magnification near the fold caustic 
is approximately A oc x~ l l 2 (Schneider, Ehlers, & Falco 1992), where x is the distance 
from the caustic, and the magnification diverges at x = (at the caustic). For a finite 
source of radius p, the magnification deviates from this form for distances x ~ 2p, and the 
finite size of the source must be considered. Binary lenses also have cusps, points where 
two fold caustics merge, and the magnification structure near a fold is quite different from 
that near a cusp. In particular, it cannot be described by the same equations. While the 
probability of encountering a fold caustic is nearly unity, the probability of encountering 
a cusp is ~ Np, where N is the number of cusps. For typical binary lenses, N = 6, and 
thus the probability is ~ 20%. If 100 events were discovered towards the bulge per year, 
and 5% of these were binaries, we would expect only one cusp crossing per year, while we 
would expect 5 events with fold crossings. Hence, we will consider only fold crossings here. 

Before continuing to the next section, where we consider the magnification one expects 
when a finite source crosses the two types of caustics considered above, we include a 
brief discussion concerning notation. Note that there are two characteristic scales in this 
problem, R e and p, which are related by equation (2.3). Either of these could be used 
as our fiducial scale. For definiteness, we choose to scale all quantities by R e . For the 
remainder of the discussion, we will also assume that t = is the time when the center of 
the source crosses (or comes closest to) the caustic. For the point-mass case, this simply 
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means setting to = 0. Finally, for the binary lens, the caustic crossings come in pairs. 
Throughout we will be referring to the second caustic crossing (since it is this crossing 
that can be predicted, see § 5), when the source is moving from the inside of the caustic 
structure to the outside. 



2.2. Extended Sources 



The magnification of an extended source is given by: 

„ , , Jd 2 rA(t;r)I(r) 

where I(r) is the intensity profile of the source, A is the magnification of a point source 
at r, and the integral is over the area of the source. The numerator and denominator are 
respectively the lensed and unlensed flux of the source. Assuming that the intensity profile 
has azimuthal symmetry, and using polar coordinates, this becomes, 

_ JP drrI(r)A(t;r) 
Aes[t) ~ J p drrl(r) ' 1 &j 

where A(r) is the angle averaged magnification function, 

A{r) = -!- j* A{r,0)dO. (2.6) 
2tt Jo 

This definition is useful because the magnification geometry of the lens is entirely isolated 
in this function. Once the function has been calculated, it can be convolved with any 
source intensity profile to give the total magnification A es . Furthermore, it is the shape of 
this function that determines how well one can resolve the source, i.e. if A is highly peaked 
at a particular value of r, then the majority of the lensed flux will be coming from a small 
range of radii near r, and therefore the lens is resolving the source. 



2.3. Angle Averaged Magnification Functions 



For both fold and point caustic crossings, the approximate form of the function A 
can be calculated analytically. When a source crosses a fold caustic, two images appear 
or disappear. For typical total binary-lens mass, the size of the source, p, is considerably 
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smaller than the Einstein ring radius of the lens, p <C 1, and thus the magnification of the 
other three images changes very little as the source crosses the caustic. Also, the curvature 
of the caustic is typically very small across the source. In this regime, the magnification in 
the vicinity of the fold caustic can be approximated by (Schneider et al. 1992), 



A(x) 



l£ 

. -tip. 



1/2 



x 1 ^ 2 , x > (inside the caustic) 



(2.7) 



where Aq is the total magnification of the three unrelated images, x is the distance 
to the caustic, and bo describes the scale of the caustic. For approximately equal mass 
binaries, bo ~ R e . For x < (outside the caustic), A(x) = Aq. Defining z = x/r, where x 
is now the distance from the center of the source to the caustic, and setting bo = R e , the 
angle-averaged magnification function is, 



A(r) = A + r- 1 / 2 j(z), 



where, 





2 l/2 
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K 



7T V 



1/2 K 



1/2 



(A) 



1/2 



Z < -1 
\z\ < 1 

z > 1 



(2.8) 



(2.9) 



Here K is the complete elliptic integral of the first kind. Figure 1 shows A(r/p) versus 
rjp for Ao = 0, p = 1, and x = —t/t e = zp = 5/3/9, p, 2/3p, 1/3/9, 0, — 1/3/?, and —2/3/9. 
For x = —t/t e > P (when the source is entirely contained within the caustic), the gradient 
of the magnification across the face of the source is small, and furthermore A(r/p) has 
no maximum. This implies that these times are not useful for resolving the source. For 
< x = —t/t e < p, A(r / p) shows a sharp peak at r = x, and thus the source is being 
resolved. However, other annuli are being significantly magnified and are thus contributing 
significantly to the total light, and therefore the resolution will be degraded somewhat. By 
contrast, for —p<x = —t/t e < 0, and in the limit — > 0, only those annuli that have 
just crossed the caustic contribute to the total light. It is therefore these times that are 
most useful for resolving the source. 

For a source of uniform brightness, the total magnification is (Schneider & Weiss 
1987), 

A es {r) = Ao + r- 1 / 2 J(z), (2.10) 

where 

z < -1 

1/21 



J(z) = 



2£/2 
3tt 



\-z)K 
&(1 + *) 1/2 (zE 



l+£ 

2 



1/2 



+ 2zE 



1/2 



(z - l)K 



(-*- 



1/2 



\z\ < 1 
z > 1 



(2.11) 
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and E is the complete elliptic integral of the second kind. Note that in equations (2.8) and 
(2.10), the dependences on z and r are separable for Aq = 0. This implies that, for a fixed 
value of z, both A and A es scale simply as p -1 / 2 . 

We now turn to the point lens case. When the separation between the lens and the 
source is much smaller than the Einstein ring radius, x <C R e , the magnification can be 
approximated by, 

A(x) = x~ 1 . (2.12) 
The angle-averaged magnification function is then, 

A{r) =r~ 1 b(z), (2.13) 



where, 



b(z) = -{l + z)- 1 K 

TT 



4z 



1 + * 



(2.14) 



and K is the complete elliptic integral of the first kind. Figure 2 shows A(r/p) versus r/p 
for p = 1, and trajectories with x = (t 2 /t 2 . + /3 2 ) 1//2 , for impact parameters (3 = (solid 
curves) and (3 = 0.5/9 (dashed curves), and t/t e = —p, —2/3/9, — 1/3/9, 0, 1/3/9, 2/3/9, and 
p. Here A achieves a local maximum whenever r = x, implying that these radii are being 
partially resolved at these times. Unfortunately, for any impact parameter (3 > 0, there is 
a range of source radii at which A is never at a maximum, those for which r < (3. These 
radii are not resolved during the caustic crossing, and thus very little information about 
the radial intensity profile will be gained for this range of source radii. 

The total magnification for a source of uniform brightness is (Schneider et al. 1992), 

A es {r) = r- l B{z), (2.15) 

where 

B{Z) = { fz [E(l/z) - (1 - z-*)K(l/z)] iVl- (2 ' 16) 

As in the fold case, the dependences on z and r in equations (2.13) and (2.15) are separable. 
Thus, for a fixed value of z, A and A es scale as p^ 1 . 



3. Error Analysis 



The flux measured from a star of finite size being lensed is F(r) = A £S Fq + B, where 
the general form of A es is given in equation (2.4). For most caustic crossing events, the 
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magnification is high, and A cs Fq >> B, thus we will henceforth use the approximation 
that B = 0. Consider a star divided up into N r radial bins, each located at r\. Defining 
Arj as the width of this bin, the flux of the star is given by, 

N r 

F(t) = ^2A{r i )I(r i )r i Ar i , (3.1) 

i 

where A is the angle averaged magnification in bin i. For the case of a large number of 
bins, or Ar^ <C p, the magnification in each bin is well approximated by equations (2.8) 
and (2.13). For a small number of bins, however, these forms are not good approximations, 
and one should use the total magnification in each bin, which can be obtained by either 
integrating equations (2.8) and (2.13) over the radii of the bin, or, equivalently, by using 
equations (2.10) and (2.15), and subtracting the flux within the inner radius of the bin 
from the flux within the outer radius of the bin, and dividing by the area of the annulus. In 
practice, the latter is simpler. The parameters one would like to recover are I(ri), the mean 
intensity in bin i for a certain wavelength. In order to do this, however, one must first know 
A(ri), which implies knowing the parameters p, (3, t e , and to f° r the point-lens case, and 
the additional parameters 6, b, and q for the binary case. In principle, one could determine 
these parameters entirely within the context of the spectrophotometric measurements. 
However, variations in the seeing and extinction make this difficult. Fortunately, this is 
not a significant hindrance because independent broad-band measurements can be used 
to constrain these parameters. Since the follow-up collaborations already monitor these 
events, and, in one case, have already used broad-band measurements to constrain these 
parameters and measure the limb-darkening of the source star (Albrow et al. 1998b), 
this should not be a problem in practice, although it does imply that the spectroscopic 
and photometric follow-up collaborations should closely coordinate their efforts. Now 
suppose that a series of spectrophotometric measurements F(tk) are made at times tk, with 
uncertainties Ofc, and these measurements are fitted to equation (3.1). The parameters of 
this fit are Ij, and the covariance matrix of the errors in these parameters is given by Cij, 
where 

c = b , b l]= ^a k — — , (3.2) 

k 1 J 

and dF{tk)/dIi = 2A(r,i)riAri. The variances in the parameters are then the diagonal 
elements of c^-, and thus the fractional error in each parameter is 8Ii/Ii = (ca) 1 / 2 / Ii. 
Assuming photon-limited precision, the fractional errors scale as the square root of the 
total number of photons received, which in turn scales as (A es Fs) 1 / 2 . From equation 
(2.10), A es oc p^ 1 , for the point caustic, and from equation (2.15), A es oc p -1 / 2 for the 
fold caustic. Let 7 be the rate at which a telescope receives photons from the unmagnified 
star in a certain spectral range, and assume that measurements are made continuously 
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during a source-radius crossing time t C ross- Then the square root of the total number of 
photons scales as /5 _iy (7t cr0 ss) 1 ^ 2 ) where v = 1/2 for the point caustic and v = 1/4 for 
the fold caustic. The fractional errors will also depend on the spatial resolution of the 
recovered intensity profile, iV r . Assuming Poisson statistics and no correlations between 

1/2 

the parameters Ii, one would expect the fractional errors in Ii to scale simply as N r ' . 
However, we find that the parameters Ii are correlated, and the errors in fact scale 

3/2 

approximately as N r . We justify this assertion in the next section. Taking these scalings 
out of Cjj, we have 

S Jl = N 3/2 htcrossr l/2 p v {dhi) l/^ (3 .3) 

where (c^) 1 / 2 is the normalized fractional error for N r = 1, p = 1, 7t C ross = 1- The specific 
form of (£j i) 1 / 2 will depend on the details of the caustic encounter and will therefore vary 
on an event-by-event basis. In the next section, we consider specific cases, and evaluate 
(cu) 1/2 directly. 



4. Results 
4.1. Scaling with Number of Radial Bins 

3/2 

In order to justify the assertion that the errors scale as N r , we consider a specific 
example for each caustic case. For the point case, we consider an event with impact 
parameter (3 = 0, and we assume that this event is observed from t/t e = to t/t e = p, 
i.e. for one source crossing time. We then evaluate the normalized error, (c^i) 1 / 2 for this 
event as a function j-j, the radius of the bin, for different values of N r . This is shown in 
Figure 3a, for N r = 10,20,40. If the errors scaled exactly as iVr^ 2 , these curves would be 
indistinguishable. In fact, the normalized errors decrease slightly for larger iV r , indicating 

3/2 

that the errors do not increase quite as fast as iV r , although this depends on rj. For the 
fold case, we assume that the event is observed from x = — t/t e = when the center crosses 
the caustic until x = —t/t e = —p when the source moves completely outside the caustic. 
In Figure 3b, we show (cu) 1 ^ 2 for N r = 10, 20,40. Here the scaling is nearly perfect. We 
will henceforth assume that the scaling is perfect for both cases, and therefore that cu has 
no dependence on iV r . 



4.2. Fractional Errors 
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In order to evaluate the fractional errors explicitly, we adopt some fiducial parameters. 
For a giant in the bulge, and a typical microlensing event, p ~ 0.02, and we assume that 
observations are taken continuously for a time, t cr0 ss ~ 7 hours. We use iV r = 10 for our 
fiducial number of radial bins, and, for the moment, we assume a uniformly bright source, 
as this is the most general way to quantify the error in a model-independent way. In 
the next section, we consider a source with limb darkening. The fractional error in the 
intensity profile for these parameters is, 

*R = (*R\ ( P V /Across V 1/2 / 7 V 1/2 [N r \ 3/2 UU 

h UJo V0.02j Uhrj U-4S-V \l0j ' l ' j 

where v = 1/2 for the point caustic and v = 1/4 for the fold caustic, and 7 = 0.4s" 1 
is approximately the flux of photons per spectral resolution element from a star with 
V = 17 that can be acquired with a 2m telescope and a spectrograph with 1 A resolution. 
In Figure 4a we show {81 as a function of r» for the point-lens case, for t/t e = to 
p, and for four different impact parameters, (3/p = 0,0.25,0.50, and 0.75 (solid, dotted, 
short dashed and long dashed lines). It is clear from this diagram, that for r > (3, the 
fractional error is small, ~ 10%. One could, in principle, improve the fractional error by 
a factor 2 1 / 2 by observing from — p < t/t e < p, i.e. for two source crossings. However, 
this will prove difficult in practice (see § 5). While the errors are small for r > (3, there is 
almost no information for radii r < (3, as expected, since the lens is not passing over this 
region of the star. In Figure 4b we show (5I/I)q as a function of r{ for the fold caustic 
with Aq = 0, and four different cases. The short dashed line corresponds to measurements 
taken from x = — t/t e = p to x = 0, one crossing time from when the edge of the source 
first crosses the caustic to when the source is exactly centered on the caustic, i.e., the first 
half of the caustic crossing. The solid line is for x = to x = —t/t e = —p (the second 
half of the caustic crossing); the dotted line is for x = — t/t e = p to x = — t/t e = —p (the 
entire caustic crossing), and the long-dashed line is for x = — t/t e = 2p to x = — t/t e = —p. 
From this figure one learns three things. First, the errors are reasonable, 57/7-20% for 
typical parameters. Second, times immediately after the center of the source crosses the 
caustic are the most crucial for recovering the intensity profile accurately. Finally, it is 
only the times just before until just after the caustic crossing that are useful for resolving 
the source. We also show the errors for the case that Aq = 1 (dashed-dot line). Clearly, 
the magnification of the images not associated with the caustic does not greatly affect the 
resultant errors. 



4.3. Effects of Limb Darkening 
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In this section we examine the effect limb darkening has on the fractional error in the 
recovered intensity profile. We can anticipate that limb darkening will serve to increase the 
errors in the outer annuli relative the constant surface-brightness case because there will be 
fewer photons coming from outer annuli, and thus the errors will be larger. For the inner 
annuli, the fractional errors will be smaller compared to the constant surface-brightness 
case, because the source is, in essence, more compact, and thus the net magnification will 
be larger, and therefore the errors smaller. 

To quantify this effect, we apply the same formalism as in § 4.3, except that we adopt 
the following parameterization of the surface brightness as a function of radius, 



where n\ and K2 are the limb-darkening coefficients. We adopt the coefficients for a cool 
(4500 K) giant (log g = 1.5) of solar metallicity from Manduca, Bell, & Gustafsson (1977) 
and Manduca (1979). These are m = 0.798,0.567,0.139 and k 2 = -0.007,0.114,0.259 for 
the V, I and K bands, respectively. Since stars are less limb darkened in the infrared, the 
results for the K band will be the most similar to the uniform surface brightness case, 
while those for V will be least similar. The results for the same fiducial parameters used 
in § 4.2 (p = 0.02, t CTOSS = 7 hours, 7 = 0.4 s -1 , N r = 10) are shown in Figure 5 along with 
the uniform surface-brightness case. For the point-caustic case [panel (a)], we have used a 
trajectory with (3 = and assumed that observations are taken for < t/t e < p. For larger 
values of (3, the fractional errors will deviate more dramatically from the uniform surface 
brightness case, since for large values of 0, only the outer annuli are being effectively 
resolved, and it is these annuli that are affected most by limb darkening. For the fold 
caustic case [panel (b)], we have assumed —p < x = —t/t e < p. As expected, for both the 
point-caustic and fold-caustic cases, the fractional errors are somewhat larger than for the 
uniform source at larger radii, and somewhat smaller at smaller radii. Also, the differences 
between the uniform source and limb-darkened source decrease for longer wavelengths. In 
general, however, for both the point and fold caustics, the differences between uniform 
and limb-darkened sources is modest. For fold caustics the effect is ~ 25% over the entire 
range of radii. For point caustics the effect is small over most of the star but the error can 
more than double at the very limb of the star when observed in V band, the bluest color 
considered here. 




(4.2) 



5. Discussion 
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While both point-mass lenses and binary lenses can in principle be used to resolve the 
two-dimensional (radial + spectral) profile of a star, binary lenses are substantially easier 
to use. First, for a binary one always has warning of the second caustic crossing. When 
the source crosses the caustic the first time it is suddenly magnified by a factor p -1 / 2 ~ 7 
and hence is easily recognized. The second crossing can then be expected in several days to 
several weeks. Intensive photometric monitoring (now routinely undertaken by PLANET 
and GMAN) can then be analyzed to make a more precise prediction. From Figure 4b it is 
clear that the most useful portion of the second caustic crossing is the final ~ 70% of the 
time that the source actually straddles the caustic. The onset of this optimal period can be 
judged extremely accurately if photometric monitoring is proceeding simultaneously, and 
reasonably accurately even one day in advance. By contrast, there is no way to guarantee 
a priori that a point-mass caustic crossing will occur because one does not know the size 
of the Einstein ring projected onto the source plane beforehand, and hence one does not 
know p. Using optical photometry alone, one can "predict" a source crossing only at about 
the time it begins. Using optical/infrared photometry, it could be predicted at r ~ 1.5 p 
(Gould & Welch 1995), but this would leave only a few hours' warning. 

Second, fold caustics generically provide information about the entire radial profile of 
the star while point caustics provide information only for annuli of the star that are greater 
than the impact parameter, (3. See Figure 4. Again, it is virtually impossible to predict in 
advance for which events (3 <C p (and hence for which events one can resolve essentially the 
whole star), although once the peak occurred, these events could be recognized provided 
that the star was being monitored photometrically. In any event, of all point-caustic 
transits, only a fraction (3/p have impact parameters smaller than [3. 

Third, binary-lens events are probably more common than point-mass caustic 
crossings. The fraction of events with binary-caustic crossings has not yet been established 
empirically, but 5% appears to be a plausible estimate. The fraction of point-mass caustic 
crossings is (p), but by the argument of the previous paragraph, only about half of these 
are really useful. The mean radius of a giant is R ~ 22 R@ (Gould 1995b), about 2.2 times 
larger than the fiducial value used in equation (4.1). Thus, the fraction of events with 
useful point-mass caustic crossings is ~ 0.5(p) ~ 2%. 

On the other hand, as shown by Figure 4, over the range probed by the point lens 
(r < (3), the point-lens errors are less than half as large as those of the binary. This 
advantage diminishes as p -1 / 4 for larger stars, but is still substantial for most giants. 
The problem of recognizing events with (3 <C p sufficiently early to permit spectroscopic 
monitoring is formidable. However, if they can be monitored beginning at their peak, they 
are the best events to follow. 
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Fig. 1. — The angle averaged magnification, A as a function of the radius in units of the 
source radius r/p for the fold caustic case, for seven different times, x = —t/t e = 5/3/9 (top 
panel), p, 2/3p, 1/3/9, 0, —1/3/9, and — 2/3p (bottom panel). The associated geometry is 
shown to the right of each panel. We have set the magnification not associated with the 
caustic to zero, and have taken p = 1.0. For other source radii, A scales as p~ 1 / 2 . 
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Fig. 2. — The angle averaged magnification, A as a function of the radius in units of the 
source radius r/p for the point caustic case, for seven different times, t/t e = —p (top panel), 
—2/3/9, — l/3p, 0, l/3p, 2/3p, and p (bottom panel), and two different trajectories, j3 = 
(solid lines) and (3 = 0.5/? (dashed lines). The associated geometry is shown to the right of 
each panel. We have taken p = 1.0. For other source radii, A scales as p^ 1 . 
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Fig. 3. — (a) The normalized fractional error in the recovered intensity profile (81 /I) 
as a function of the radius in units of the source radius r/p for the point caustic. Each 
line corresponds to a different number of radial bins N r in the recovered intensity profile, 
N r = 10 (solid), 20 (dotted), and 40 (dashed), (b) Same as (a) except for the fold caustic. 
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Fig. 4. — The fractional error in the recovered intensity profile (51 /I) as a function 
of the radius in units of the source radius r/p for the fiducial parameters: N r = 10 
(number of radial bins in the recovered intensity profile), 7 = 0.4 s -1 (unlensed flux of 
star), tcross = 7 hours (crossing time of star) and p = 0.02. (a) The errors for the point 
caustic, for t/t e from to p and four different values of the impact parameter, [3/p = 
(solid), 0.25 (dotted), 0.5 (short dashed), and 0.75 (long dashed), (b) The errors for the 
fold caustic, for four different trajectories: x = —t/t e from to — p (solid), p to (short 
dashed), p to —p (dotted), and 2p to —p (long dashed). For these curves, we have set the 
magnification not associated with the caustic to zero (^4o = 0) . We also show the fractional 
error for the case A§ = 1 and x = —t/t e from to —p (dashed-dot). 
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Fig. 5. — (a) The fractional error in the recovered intensity profile (51 /I) as a function 
of the radius in units of the source radius, r/p, for the point caustic, for a uniform source 
(solid), and limb-darkened source as appropriate for the broad bands V (dotted), I (short 
dashed), and K (long dashed), (b) Same as (a) except for the fold caustic. For both (a) 
and (b) we have assumed the same fiducial parameters as the solid curves in Figure 4. 



